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In this work I would like to present some ideas on how to improve on the gauge sector 
in our lattice simulations at finite baryon density. The long standing problem, that we 
obtain an onset in thermodynamic quantities at a much smaller chemical potential than 
expected, could be related to an unphysical proliferation of flavours due to hard gluons 
close to the Brillouin edges. These hard gluons produce flavour non-conserving vertices 
to the fermion sector. They also produce excessive number of small instantons due to 
lattice dislocations. Both unphysical effects could increase the propagation in (di)-quarks 
to give the early onset in \x. Thus we will present here a modified action that avoids large 
fields close to the lattice cutoff. Some of these ideas have been tested for 577(2) and are 
being implemented for 577(3). 

1. Introduction 

There has been the long standing problem, since the introduction of a chemical potential 
fj, coupled to the baryon density |],[2[], that the onset of thermodynamic quantities happens 
in our lattice simulations @-[5| at a much smaller fi than at the expected \x = The 
baryon density and the chiral condensate start changing from their vacuum values at \x of 
the order of or earlier than m 7r /2, as if controlled by a Goldstone mode with net Baryon 
number, which is clearly unphysical in the low Temperature normal phase. 

The community had the hope that including properly the fermion determinant, which 
for non-zero \i is complex and therefore very difficult to implement with lattice methods, 
will by itself solve this problem. The careful work by the Barbour group in ref. || and 
the follow up presented by I. Barbour and S. Morrison at this Workshop, seem to show a 
pinching at a higher \x but the onset still seems to be around m n /2. 

Several years ago I described another possible problem, on top of the fermion deter- 
minant one, that could be producing the early onset in our lattice simulations. As is 
well known, for Kogut Susskind fermions the flavour number is not conserved along their 
propagation due to hard gluons that can make jump the fermion propagator pole to other 
corners of the Brillouin zone, thus changing its flavour. This flavour changing would 
produce a proliferation of flavours even for the valence quarks, forming a much tighter 
bound nuclear matter than in nature . As we heard on various talks by F. Wilczeck, 
E. Shuryak and others at this workshop, diquark condensates could be produced easier 
the more flavour species, which could drive the observed early onset due to the flavour 
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proliferation. Furthermore, these hard gluons, with momenta p « it /2a or higher, pro- 
duce easily lattice dislocations that give many unphysical "small instantons" || . Again E. 
Shuryak in his model of quark propagating among instantons, gets increased propagation 
with a higher instanton density. 

We see then, that from several points of view it seems important to find a way to 
suppress the gluons with momenta of the order of the lattice cutoff. We could regain 
flavour conservation and suppress unphysical mini instantons. Both unwanted effects 
tend to lower the \x onset value from the physical one. 

In order to achieve such a hard gluon suppression I have been experimenting, first for 
SU(2), with an "improved" gauge action which has the same naive continuum limit as 
the usual one in the field strengths, but diverges for plaquette phases of the order of the 
lattice cutoff. Even though this suppression method is gauge invariant, it is important 
to use a Metropolis updating scheme with small changes in the link phases, in order to 
guarantee that we stay locally without spurious dislocations. In this way we can hope 
to suppress small unphysical instantons. The physical gluon momentum transmitted to 
the fermions, given in the continuum by the Poynting vector, being essentially in the 
nonabelian case shaped as a "chair" formed with two contiguous plaquettes, could then 
also be kept low so as to conserve flavour. 

With L. Polley we had also tried out a new kind of lattice fermions P,|lOl , the t-asymmetric 



ones, where one can reduce to 2 flavours from the start (in contrast to the 4 staggered 
ones), but have several technical problems to surmount, related with regaining the 0(4) 
symmetry, in order to calculate the hadron masses and compare them with the onset. 
This will be treated in the next section. 

In sec. 3 we will describe the new gauge action, mainly for 577(2), and discuss its 
implementation for SU(3). The ultimative test will be to generate configurations gener- 
ated with such an improved gauge action, to then compute the baryon density or chiral 
condensate for the full theory. 

2. t-asymmetric fermions 

An earlier attempt in order to reduce the number of flavours, was based on a new kind 
of fermions, also called t-asymmetric fermions, that don't allow more than 2 flavours in 
contrast to the 4 flavours for staggered ones PJTO|]. This fermions are obtained by just 
taking a one-sided time derivative, which still pro duces an hermitian hamiltonian and 
eventually a positive transfer matrix (for a T < a s j ^3 + rri£). Consequently, the number 
of flavours gets reduced to just 2. These fermions are related to Susskind's hamiltonian 



formalism ones, via a unitary transformation [|il| , |12| and a related form is also being 
investigated by the Heidelberg group [13| . 



The details of this formulation can be found in ||1G], with the result that after a 



diagonalization in spinor space like the usual Kawamoto- Smit one, 

^{x) = a^afa^xix) with a k = i^^ k (1) 
we get, 



k 



3 



with the usual: T u (x) = (— ) 9! i+ , "+ a! »'-'- ) after having thinned out to the first diagonal 
component. This fermionic action xKx nas then just 2 flavours. Susskind showed that 
there is a discrete version of chiral invariance left over, preventing mass counterterms 
in the interacting case, but not guaranteeing "Goldstone behaviour" for the pions. One 
could still get light pions, like in Aoki's proposal for Wilson fermions fll4" |. 

We had already done first simulations with the baryon density for these fermions, 
(J a) = (e^Xx U x> iX x +i) an d the chiral condensate, (%x)> as function of \x. The operators 
have been calculated with the solved pseudofermion method ||]. The curves for (J) are 
shown in ref. |I| (Fig. 1) and show a clear onset for two m q = .01, .04 at the same 
\i pa 0.2. This result could be deceiving though as it is not clear how m n scales with m q 
for these fermions, as already mentioned above. In order to really check if we are getting 
better results, we have to compute for these fermions the ir and N masses. For this, 
considering that we need a finer lattice in the r direction for positivity and a fine a T for 
charge conjugation symmetry in the propagators, we have to scale the couplings to regain 
the 0(4) symmetry. Once this is done, one can attempt to calculate the masses and the 
baryon density for this 0(4) symmetric action. This has proven to be a formidable task 
and presently I consider more promissory to improve on the gauge sector. 



3. Improvements on the gauge action 

The other way out, in order to prevent flavour changing is to suppress the appearance 
of hard gluons with momenta p > ii/2a. Without these hard gluons the quarks cannot 
change their flavour by jumping to other corners of the Brillouin zone. As mentioned, 
these hard gluons can also cause unphysical mini-instantons that could enhance the quark 
propagation in Shuryak's scheme, also affecting the \x onset. 

From my investigations in SU(2) with instantons, in order to try to suppress mini- 
instantons produced by lattice dislocations, one can do so by modifying the action for 
plaquettes with a large phase. Instead of the usual gauge action: 

S = (! - l / N Re tr([/ a )), (3) 
□ 

which has a maximum plaquette action of 2f3 for a phase n = n, one can take an improved 
action of the form: 

S = j3 ]T 2/tt tan( vr/2 (1 - l/N Re ti(U a ))), (4) 
□ 

which has the same naive continuum limit for small a, but diverges for plaquettes with 
phases close to the cuttof. One can easily device functions which start the same and grow 
even faster for larger plaquette phases. The point being that at relatively moderate /3's 
one can suppress plaquette phases larger than n ~ 7r/2 and therefore make sure that the 
effective tr(F^) stays away from the cutoff. ( For the usual gauge field action in SU(3), 
for example at (3 = 6.0, around 1/3 of the plaquettes still have a phase On > tt/2 !). The 
actual gauge invariant definition for the gluon momentum, resembling the Poynting vector 
(trE A B) in the continuum, involves pairs of orthogonal plaquettes like a "chair" . If each 
of the plaquette phases is kept small also the combined one should be constrainable to 
get p < 7r/2a. 
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In order to avoid large phases in the gauge links, which produce dislocations earring 
a topological charge 0, it is convenient to use the Metropolis algorithm for the link up- 
date, with small changes in the SU(N) matrix on each iteration. In combination with an 
improved action as introduced above, one can obtain then fairly "smooth" gauge configu- 
rations, starting from the identity one. There is a problem though, related with a random 
walk like drifting of the links at a vertex corresponding to local gauge transformations, 
until we get large link phases close to tt that can produce dislocations. This can be re- 
duced by gauge transformations that try to keep all the link phases low. The trick I used 
in SU(2) is to check the phase <fi of the new link U = exp(i <fi cf) ■ a) and if it is larger 
than some phase, lets say 7r/6, make a local gauge transformation that shrinks that link's 
phase to it/8 keeping its SU(2) direction, while applying the inverse element to the other 
links at that vertex. These local gauge transformations succeed in keeping almost all <j) 
smaller than tt/6, thus suppressing almost all dislocations. 

In this way I have been able to constraint the total topological charge to values 
Qv < | ± 1 1 in periodic lattices (the field-theoretic definition does not give exactly in- 
tegers), for fairly large 577(2) lattices (12 3 x 24) with f3 = 2.5. For these lattices also the 
charge Q v with v — V/2 half the volume is also smaller than | ±1| showing no I — I pairs. 

One could also use more conventional action improvements in order to avoid hard 
gluons, by including terms in the action proportional to the Poynting operator (more 
generally the energy-momentum tensor 0^) that suppress large fields. In other words, 
include "chair" diagram terms in the action. 

Another interesting alternative for an improved action in order to reduce the flavour 
breaking, was presented by Lagae and Sinclair ||15| . They smear the links by covariant 



displacements in all directions, in order to project the momenta to the central pole region 
in the Brillouin zone. In momentum space in the continuum they project the fields 
Afj,(k) — > 1/16 TI i=1 4 (1 + cosfcj) A M (fc), which in coordinate space corresponds to the 
smearing A^x) — > 1/256 Ili=i..4(2 + A: + D_i) A^(x). They report a large reduction 
in the flavour symmetry violations for the (non)-Goldstone ir splitting of masses. Up to 
now the only introduce this projected gauge fields in the minimal coupling terms in their 
action and it would be interesting to see what happens with the topology if they consider 
such improvement in the gauge action. 
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4. Conclusions 

The inclusion of the fermion determinant could be not the only ingredient needed in 
order to shift the onset chemical potential to physical values close to the Nucleon mass. 
Our lattice simulations at finite \x have been done with an unimproved gauge sector. This 
produces several unwanted effects, like the proliferation of flavours even for valence quarks 
due to the hard gluons that can change the flavour along the quark propagation. Also 
the topological sector is distorted due to unphysical mini-instantons appearing due to non 
smooth gauge fields. 

Both unwanted effects should be suppressed with an improvement in the gauge action. 
I introduced a simple to implement improvement that should help avoid hard gluons and 
dislocations. These ideas have been tested for SU (2) with encouraging results and could 
also be implemented for 577(3). In the Metropolis algorithm we cannot use anymore the 
trick to have the sum of staples precalculated as the action is not linear anymore, but the 
rest is very simple to implement. 

Once we have the gauge sector under control, with no hard gluons close to the cutoff and 
no unphysical mini-instantons, we could again calculate the baryon density operator or the 
chiral condensate to find its /x dependence. This gauge updating could also be integrated 
in the unquenched code to redo Barbour's group method to include the determinant. 

I would like to thank L. Polley for several interesting conversations on this gauge field 
improvement. 
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